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Several sources, both experimental and theoretical, subject inclusive jet cross section measure- 
ments to large uncertainties. On the experimental side, the energy scale contributes between 10% 
and 30% uncertainties. On the theoretical side, choice of parton distribution function introduces a 
20% variation in the prediction; freedom in renormalization scale provides another uncertainty of up 
to 30%. The ratio of inclusive jet cross sections reduces the uncertainty from these major sources, 
permitting a very precise test of next-to-leading order QCD. The preliminary results from D0 and 
CDF differ from simple QCD predictions in normalization and differ from each other at small values 
of jet xt- 
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At sufficiently high energy, field theory calculations reduce the complexities of proton-antiproton interactions to 
simple scattering processes involving only one constituent from each particle. Identifying these constituents, called 

I ^~ , "partons," as quarks and gluons, quantum chromodynamics (QCD) calculations determine the matrix elements for 

,-4- | a given scattering process, but require input from empirically-determined parton momentum distribution functions 
(pdf's). Perturbative QCD calculations provide specific predictions for the rates of scattering processes as a function 
of final-state parton momentum and angle, but the unknown initial-state momentum superimposes a boost to the 

Q\ , final system. Consequently, cross section predictions reflect an average over the initial momenta through sampling 
the pdf over a large number of events. The two largest uncertainties in the QCD prediction result from the choice of 

iq ■ pdf and the choice of renormalization scale. 

Experimentally, a final state parton is observed as a shower of collimated particles: a jet. Jet cross sections 



are reported in terms of Et, the jet energy in the plane transverse to the beamline. This analysis compares the 
production rate of jets at two center-of-mass energies, 630 and 1800 GeV, for two reasons. First, the ratio of cross 



sections constrains experimental uncertainties that are common to both studies. Second, the ratio suppresses the 



QCD prediction's sensitivity to the choice of pdf. The ratio of cross sections thus provides a stronger test of the 
QJ \ matrix element portion of the calculation than a single cross section measurement alone. 
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•i-l , II. CONTEXT 
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A simple parton model without gluon bremsstrahlung would predict a jet cross section that scales with center-of- 
mass energy (y/s). In this scenario, a dimensionless jet cross section, Ej, ■ Ejr-%-, plotted as a function of jet xt = ^7=, 
would remain constant with respect to changes in center-of-mass energy. 

Although past cross section ratio analyses at lower y/s exhibited significant deviation from the scaling model, the 
dimensionless framework continues to provide a useful context for comparisons with the more developed QCD theory; 
thus, the prior tests motivate the selection of variables in this work. The time-integrated luminosity provides a 
measure of the statistical power of these earlier cross section measurements: the 1985 results from UA2 pi include 
122 nb" 1 (y/s = 546 GeV) and 310 nb" 1 (y/s = 630 GeV) of data, whereas the 1993 result from CDF § compared 
8.58 nb- 1 (y/s = 546 GeV) to 4430 nb" 1 (y/s = 1800 GeV). 

As reported in jj|, the D0 collaboration at Fermilab measured the inclusive jet cross section at y/s — 1800 GeV using 
93,000 nb -1 of data. This article presents a complementary inclusive jet cross section analysis at y/s — 630 GeV, with 
537 nb -1 of data. The most recent analyses from the CDF collaboration at Fermilab use data samples of similar size 
at both center-of-mass energies. Because in each case, the data sets for both center-of-mass energies were collected 
with the same detector B B, many uncertainties in the results are highly-correlated, particularly uncertainties in the 



calorimeter response and calorimeter noise. Additionally, the explicit expression of the correlations in uncertainties 
from bin to bin provides a more stringent comparison to QCD than was previously possible. 

The remainder of this article focuses on the methods of the D0 analysis of the jet cross section and ratio, reserving 
discussion of the CDF result for last. Although the details differ, both experiments require the same components for 
their respective analyses: data reconstruction, sample selection, energy scale correction, and resolution correction. 

III. DATA RECONSTRUCTION AND SAMPLE SELECTION 

The D0 collaboration employs the variables transverse energy (Et) and pseudorapidity (n = — In (tan |)) to de- 
scribe the hadronic particle showers that result from inelastic pp collisions. Here, 9 measures the polar angle relative 
to the proton beam direction. Jets are defined by their hadronic shower Et and r\ — (j) centroid. Starting with the 
most energetic deposits in the calorimeter, the D0 reconstruction algorithm U iteratively finds the energy-weighted 
centroid for each hadronic shower within a cone of dimensionless radius 0.7. When two such showers overlap, they 
are merged into a single jet if they share more than 50% of the Et of the lowest-£/r shower; otherwise, they are split 
into two separate jets, each with its own centroid and Et- 

The online data selection procedure triggers on events that contain at least one jet above threshold. A small 
correction, 3% for the first data point and decreasing thereafter, removes the effect of a slight trigger inefficiency at 
low Et- The offline data selection procedure, which eliminates background caused by electrons, photons, noise, or 
cosmic rays, closely follows the methods J7| || of the 1800 GeV analysis. The efficiency of jet selection is nearly 
constant as a function of jet Et, approximately 96%. To maintain optimal precision in jet energy measurements, a 
vertex requirement discarded jets originating from interaction points more than 50 cm from the longitudinal center 
of the detector, lowering the total data selection efficiency to 82%. The cross section uncertainty associated with all 
efficiencies amounts to less than 0.5%. 

IV. ENERGY SCALE AND RESOLUTION CORRECTIONS 

The jet energy scale correction, described in 0, corrects the energy response of the D0 calorimeter and removes 
cone boundary effects, multiple interaction effects, calorimeter noise, and the underlying event energy resulting from 
spectator partons. The response correction increases the Et of jets by 22% when the raw calorimeter Et is 20 GeV, 
and increases the Et by 10% at asymptotically-high raw Et's. Calorimeter noise contributes on average 1.6 GeV 
of Et to each jet. The underlying event, which represents the only significant difference between center-of-mass 
energies, contributes 0.6 GeV to each jet at y/s = 630 GeV, compared to 0.9 GeV at y/s = 1800 GeV. Simply 
speaking, a jet's measured Et increases by 12 to 14% after the total energy scale correction. Uncertainties from the 
noise and response corrections dominate the systematic uncertainty band of the final result. 

Although the energy scale algorithm corrects jets from their reconstructed Et to their "true" Et on average, 
the energy scale cannot remove the random fluctuations of individual jets about this average value. The resulting 
imperfect resolution of a jet's Et about its true value leads to a smearing effect that effectively inflates the observed 
cross section, especially in the steepest portions of the distribution. After removing contributions that do not result 
from detector effects, deviations from perfect momentum balance in dijet events provide a measure of the average 
resolution as a function of jet Et- 

The observed cross section distribution can be described as the convolution of the resolution function and an initial 
distribution, F(Et) — f G(Et — E' T ) ■ f(E' T ) dE T , where G is a Gaussian distribution and the initial distribution / 
is modelled by: 

f(E T )=A(E T ) B (l-^) D . (1) 



The Gaussian has -E/r-dependent width, and represents the jet Et resolution. The three parameters of the ansatz 
function, A, B, and D, are determined with a best-fit of the smeared distribution, F, to the observed cross section. 
The correction factor for the smearing effect is expressed as C(Et) — J{Et)/F{Et); the measured cross section is 
multiplied by C on a bin-by-bin basis. The uncertainty of unsmearing the cross section is expressed in an error matrix, 
where the partial derivatives of the correction factor with respect to each parameter are determined numerically. In 
addition to the fitting uncertainties of C, this matrix includes contributions from the resolution determination. 
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FIG. 1. (Left) The inclusive jet cross section at ,/s = 630 GeV. Black points indicate the cross section measurement and 
statistical errors; histogram line indicates NLO QCD Prediction from Jetrad. (Right) Fractional deviation of the cross section 
data and NLO QCD. Shaded blocks indicate systematic uncertainty. 



V. RESULTS 



The unsmeared y/s = 630 GeV jet cross section appears in the left side of Figure [l| The point position within each 
Et bin occurs where the ansatz takes its average value within the bin. The corresponding xt positions appear on the 
upper axis in anticipation of the ratio of cross sections; by design, the xt points provide an exact match to the points 
derived from the jet cross section at y/s = 1800 GeV. The right side of Figure [j] compares the cross section to a NLO 
QCD prediction, with renormalization and factorization scales set to Et/2 and using CTEQ3M parton distribution 
function. R se p, a phenomenological parameter jlQ] in the prediction, was set to the usual value of 1.3. The shaded 
blocks in the figure indicate the la systematic uncertainty of the cross section measurement, while the vertical bars 
indicate the la statistical uncertainty. 

The primary uncertainty in the ratio of cross sections results from the energy scale correction, and the remaining 
uncertainty consists mostly of the resolution and luminosity uncertainties, as depicted in Fig. 0. Although the 
systematic error in the individual cross sections ranges from 10% to as much as 30%, strong numerator-to-denominator 
correlations in each bin of the ratio of cross sections reduce the uncertainty to as low as ±5%. 

Because they result primarily from best-fits, the uncertainties associated with the energy scale and unsmearing 
corrections do not allow the cross section ratio to normalize freely. Instead, deviations in one portion of the distribution 
increase the probability for deviations in other parts of the distribution. The correlation of one bin in the ratio to 
another determines the shape evolution of the ratio, given a proposed displacement at a particular point. The 
correlation information is represented in the bottom of Fig. 0, where the diagonal elements are unity by definition. 
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FIG. 2. (Top) Fractional uncertainty in the ratio of cross sections, displayed by component as functions of jet xt- (Bottom) 
Correlations in the uncertainty from bin to bin. 

The observed ratio of scaled cross sections takes a value between 1.4 and 2.0, depending on the position in xt- As 
shown in Fig. 0, NLO QCD predictions for the ratio of cross sections lie systematically above the data throughout the 
majority of the measured xt range, most notably in the range near 0.17, where the ratio has the most statistical power. 
Many theoretical choices, including choice of parton distribution function, minutely change the ratio result; only the 
renormalization/factorization scale changes the prediction significantly. The right-hand side of Fig. |3| compares the 
data to predictions with various renormalization scales, which were generated with Jetrad pl|. 
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FIG. 3. The ratio of dimensionless jet cross sections, with NLO QCD predictions for: (left) various pdf 's and (right) various 
renormalization scales. 
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FIG. 4. The ratio of dimensionless jet cross sections, versus predictions with different [i, scales in the numerator and the 
denominator. 



PDF 


fi Scale 630 GeV 


fi Scale 1800 GeV 


x 2 


prob. 


CTEQ3M 
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jEt 
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Et 
^Et 
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Et 
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Et/2 
Et/2 
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10.9 
28.5 
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9.84% 
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TABLE I. D0's preliminary \ comparisons for the ratio of jet cross sections (20 degrees of freedom). 



The matrix of errors and correlations allows a x 2 comparison between data and theory that can quantitatively 
describe the probability that the two curves describe the same parent distribution. The number of data points, 20, 
determines the number of degrees of freedom. The results, in Table M, indicate poor agreement between the data 
and NLO QCD for predictions with the same renormalization scale (/i scale) in the numerator and the denominator. 
For predictions with different /j, scales for two y/s, agreement improves to as much as a 95% probability that the 
data and prediction represent the same distribution (Fig. 0). The need for different \i scales may indicate a residual 
dependence of the prediction to higher-order terms in the calculation. 

VI. THE PRELIMINARY CDF RESULT 

Shortly after DPF '99, the CDF collaboration updated its preliminary result with two major improvements. First, 
the ratio of cross sections reflects an improvement in the available time-integrated luminosity of the 1800 GeV cross 
section, which increased from 4,430 nb _1 to 87,000 nb _1 . (At 630 GeV, the updated CDF luminosity totals 576 nb _1 .) 
Second, the result |1^] now includes a systematic error estimate (Figure g) . 

A comparison of Figure y| and Figure |5j reveals excellent agreement between the experiments above jet xt of 0.12, 
but sharp disagreement below that value. Currently, the source of this disagreement between experiments is not 
understood. The CDF data in Figure [| are entirely consistent with the CDF result of 1993 in Reference [||. 
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FIG. 5. CDF's preliminary ratio of dimensionless jet cross sections. 



VII. CONCLUSIONS 



The ratio of inclusive jet cross sections provides a strong test of NLO QCD matrix elements while suppressing the 
uncertainties that result from pdf's. These D0 and CDF results are drawn from larger data sets than previously 
available; thus, they span a larger range of jet xt and exhibit smaller uncertainties than any prior measurement. No 
quantitative comparison is needed to see that the present CDF result exhibits a clear disagreement with D0 at low 
values of jet xt- 



The x 2 tests of D0 indicate that NLO QCD does not model the D0 data unless a different renormalization scale 
is introduced for the two center-of-mass energies. Because the renormalization scale has no physical meaning, the 
broad range of % 2 probabilities suggest that the precision of the data exceeds the precision of the NLO predictions. 
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